Angle dependent magnetoresistance measurements in Tl2Ba2Cu06+5 and the need for 

anisotropic scattering 
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The angle-dependent interlayer magnetoresistance of overdoped Tl2Ba2Cu06+i has been mea- 
sured in high magnetic fields up to 45 Tesla. A conventional Boltzmann transport analysis with no 
basal-plane anisotropy in the cyclotron frequency uJc or transport lifetime r is shown to be inade- 
quate for explaining the data. We describe in detail how the analysis can be modified to incorporate 
in-plane anisotropy in these two key quantities and extract the degree of anisotropy for each by 
assuming a simple four-fold symmetry. While anisotropy in ujc and other Fermi surface parameters 
may improve the fit, we demonstrate that the most important anisotropy is that in the transport 
lifetime, thus confirming its role in the physics of overdoped superconducting cuprates. 

PACS numbers: 74.72.Jt, 71.18.-|-y, 72.10.Bg, 73.43.Qt 
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I. INTRODUCTION 

There are many routes to investigating the mechanism 
of high-temperature superconductivity and naively one 
might expect the normal state to be the simplest. De- 
spite concerted experimental effort however,^ the nor- 
mal state properties of cuprates remain a profound the- 
oretical challenge.^ Indeed, even from the earliest trans- 
port measurements in these compounds it was clear that 
the normal state was far from conventional."^ Arguably 
the most remarkable phenomena are the distinct power 
laws of the in-plane resistivity pab and inverse Hall an- 
gle cot0H temperature dependences. In optimally doped 
YBaaCuaOr-s (YBCO) and La2-:.Sr^Cu04 (LSCO) for 
example, Pab{T) varies linearly with temperature over 
a wide temperature range, whereas cotOu maintains a 
strong dependence. ^'^ In other words, it is as if these 
materials exhibit distinct scattering mechanisms which 
are separately manifested according to the experimental 
probe being considered. Anderson coined the phrase 'life- 
time separation' to describe this anomalous behavior and 
today its interpretation remains one of the greatest obsta- 
cles to the development of a coherent description of the 
normal state quasiparticle dynamics in high-Tc cuprates. 

Three contrasting approaches dominate the current 
thinking on the transport problem in cuprates; Ander- 
son's two-lifetime picture,^ marginal Fermi-liquid (MFL) 
phenomenology^ and models based on fermionic quasi- 
particles that invoke specific (anisotropic) scattering 
mechanisms within the basal plane. ^""'^^ In the two- 
lifetime approach, scattering processes involving momen- 
tum transfer perpendicular and parallel to the Fermi 
surface are governed by independent transport and 
Hall scattering rates l/rtr and I/th with different T- 
dependencies. The proponents of the MFL hypothesis 
assume a single T-linear scattering rate which naturally 
accounts for pab{T), but introduce an unconventional ex- 
pansion in the magnetotransport response whereby the 
Hall angle, for example, is given by the square of the 
transport lifetime. This anomalous expansion is at- 



tributed to anisotropy in the (elastic) impurity scattering 
rate, possibly due to small-angle scattering off impurities 
located away from the Cu02 plane. 

Attempts to explain the anomalous behavior of pab{T) 
and cot0H(T') in cuprates within a Fermi-liquid (FL) ap- 
proach have centered around the assumption of a (single) 
inelastic scattering rate that is strongly dependent on 
the quasiparticle wave-number k. This anisotropy can 
arise either due to anisotropic electron-electron (possibly 
Umklapp) scattering^^ or coupling to a singular bosonic 
mode, be that of spin,^'^ charge"'^" or d-wave supercon- 
ducting fluctuations.^^ Generating a clear separation of 
lifetimes within these single-lifetime scenarios however re- 
quires a subtle balancing act between different regions in 
k-space with distinct T-dependencies.^* 

In order to test these various proposals and to pro- 
ceed towards a theoretical consensus, information on the 
momentum (k) and energy (w or T) dependence of the 
transport lifetime r at or near the Fermi level ep is 
urgently required. This is a non-trivial exercise how- 
ever since the transport coefficients themselves are angle- 
averaged quantities involving differently weighted inte- 
grations around the Fermi surface (FS). Whilst angle re- 
solved photoemission spectroscopy (ARPES) can probe 
directly the in-plane quasi-particlc lifetime via the imag- 
inary part of the self-energy ImE(k,tt'), its relevance to 
DC transport is still unclear. ^'"^ Moreover, there remains 
some dispute as to the correct form of ImE(k,ti;) even for 
samples with nominally the same composition. ^^^7 

Measurements of interlayer magnetoresistance as a 
function of angle have yielded important information 
about the FS topology (size and shape) in a variety of 
layered metals including organic conductors'^ and quasi- 
two-dimensional (q2D) oxides. '^"^' In a recent paper, we 
showed that this technique could be developed to extract 
information on the scattering rate anisotropy and applied 
the technique to overdoped TlaBaaCuOe+a (T12201).22 
In the present paper, we present a more thorough and de- 
tailed analysis of our angle-dependent magnetoresistance 
(ADMR) measurements on overdoped T12201, focusing 
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in particular on the procedure used to fit ADMR and 
show how this analysis fails at higher temperatures un- 
less one includes such anisotropy in r. We progressively 
introduce anisotropy into the formalism, and explore the 
effects of this both in the cyclotron frequency lOc and the 
transport lifetime r. The approach presented here is sim- 
ilar to that described recently by Kennett and McKenzie 
who derived a generalized expression for ADMR in lay- 
ered metals with basal-plane anisotropy. ^"^ In this paper, 
we focus on issues pertinent to T12201, the importance 
of each parameter in fitting the ADMR signal and their 
interdependence, and the issue of sample misalignment. 
Although it is difhcult to isolate anisotropy in one from 
anisotropy in the other, the strong temperature evolu- 
tion of the ADMR signal (and subsequent measurements 
of its doping dependence^^) suggests that the dominant 
anisotropy is in r and not ujc- The paper is set out as 
follows. Section II describes the FS parameterization of 
T12201 and the necessary symmetry considerations with 
respect to the ADMR analysis. Section III briefly de- 
scribes the ADMR experiment itself. The Boltzmann 
formalism and the resulting analysis is described in Sec- 
tion IV for the cases where the parameters Uc and r are 
both isotropic and anisotropic (within the basal-plane). 
Our conclusions are presented in Section V. 



II. THE THREE-DIMENSIONAL FERMI 
SURFACE OF TL2BA2CUO6+5 

The present FS parameterization is identical to that 
used previously^"^^^'^'^ and so shall be described here 
only briefly. The interested reader is referred to [Ref.^^] 
which details a similar parameterization of the a sheet 
of Sr2Ru04. Being extended in the kz direction, the 
quasiparticle dispersion contains a finite (though small ^ 
meV) transfer integral t±{(f>, fc^), where (f> is the azimuthal 
angle in the /cj. — ky plane. The parameter t±{4>,kz) is 
anisotropic in the plane and the Fermi wave-vector kp 
is therefore modulated by both the in-plane dispersion 
and by t±{4>, kz)- The clearest way to express this is by 
expanding kp into cylindrical harmonics^®'^° 
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FIG. 1: (Color online) The Brillouin zone stacking for 
Tl2Ba2Cu06-|-a. There is inversion symmetry about the 
kx — ky plane: n,ij> ^ —k, and three further symmetries 
about the axis along the X-hne. As described in the text 
they are (i) the two fold symmetry (j) ^ (p + tt, (ii) the mirror 
symmetry (f) — 7r/2 — (p and finally (iii) the screw symme- 
try involving a translation in the kz direction and a rotation: 

K, (p ^ K + TT, (f> + Tv/2. 
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where {g} denote coefficients corresponding to cosine 
and sine terms, k = kzc/2 and c = 23. 2A is the inter- 
layer dimension of the unit cell. The symmetry of the 
Brillouin zone limits the number of parameters of inter- 
est and a pictorial illustration of this is shown in Figure 1 . 
In the kz direction, inversion symmetry k — > — k requires 
that the only terms containing k are cosines. Three fur- 
ther symmetries restrict the parameterization: (i) the 
two-fold rotational symmetry — > -I- tt, (ii) the mirror 
plane (j) — > — 7r/2 — and (iii) the screw symmetry 



K + TT, (f) + tt/2. The transformations differ from those in 
Ref.^^ because of a different choice in coordinate axes, 
but the operations are identical. The first symmetry re- 
quires that all TO be even. The next symmetry requires 
that all cosine terms have m mod4 = and all sine terms 
have TO mod4 = 2. For example, cosAcp = cos4(— 7r/2 — cp) 
whereas cos20 = — cos2(— 7r/2 — (f>). The reverse is true 
for the sine terms. The final symmetry requires that all 
of the cosine terms be accompanied by n that are even 
and the sine terms be accompanied by any n that are 
odd. For example, cosKsin2(/) = cos(k -|- 7r)sin2(7r/2 -I- (f>), 
but cos2Ksin20 = — cos2(k + 7r)sin2(7r/2 -I- cj)). The con- 
verse is of course true for the cosine terms that have 
TOmod4 = 0. Eq. 1 can thus be simplified to 





FIG. 2: (Color online) Quasi-2D Fermi surfaces described by 
Eq. 2 with progressive inclusion of cylindrical harmonics (a) 
fcoo, (b) k2i, (c) ko4 and (d) fcei and fcioi. 
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FIG. 3: (Color online) The projection of the 3D dispersion 
for overdoped T12201 as determined by ADMR (thin black 
lines), plotted over the ARPES (thick green line) results for 
a compound with a nominally similar doping. The warping 
of the transfer integral is exaggerated by approximately 100 
times for clarity. 



III. EXPERIMENTAL 
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We have shown previously that the minimum number 
of parameters required to fit the data that simultaneously 
satisfy these symmetry constraints are ^007^047^1:^61 
and fcioi.^" Figure 2 shows the warping created by pro- 
gressive inclusion of these parameters, beginning with a 
dispersionless {t± = 0) isotropic FS. Eq. 2 has exact four- 
fold symmetry though the modulation t±_ gives rise to 
eight highly symmetric points where the transfer integral 
vanishes, as predicted by band-structure calculations.^"* 
Figure 3 shows the projection of the three-dimensional 
FS as deduced by ADMR^" overlaid on that determined 
by ARPES^^ (see Eq. 15; this curve corresponds to the 
nominal doping level of this crystal). The agreement is 
very good, but most importantly the two experiments, 
to a good approximation, share the eight points of high 
symmetry. For ease of computation, the Fermi surface 
can be described by 

eF(t,^) = ^4^(0) - 2txcosK 

x-(fc2isin20 -I- /cgismG^ -I- /cioisinlO^), 

(3) 

where fcp((/)) = fcoo + fco4cos(40) and a = 3. 866 A is the 
in-plane lattice parameter. 



T12201 is the most suitable cuprate system for ADMR 
studies due to its single Fermi sheet, its strong two- 
dimensionality,^^ its low residual resistivity^^"'^'' and ac- 
cessibility to the whole overdoped region of the cuprate 
phase diagram. "^^ Single crystals were fabricated using a 
self-flux method in alumina crucibles. As-grown crys- 
tals are naturally overdoped and the doping level (and 
therefore the desired Tc) is set by annealing in oxygen, 
argon or in vacuum. '^^ The crystal used in this study 
(300/im x 150/im x 20/im) was annealed in oxygen at 
600K for 200min resulting in « 17K. Electrical con- 
tacts were attached using Dupont 6838 silver paste in 
a quasi-Montgomery 4-wire configuration. ADMR mea- 
surements were performed at 45T in the hybrid magnet 
at the National High Magnetic Field Laboratory, Talla- 
hassee, Florida using a probe with a two-axis rotator. 
Initially, the platform on which the sample was mounted 



was rotated by an azimuthal angle 



and then the 



interplane resistivity pzz was measured as the polar an- 
gle 9cxp was swept at constant temperature and constant 
field (see Figure 4). 



IV. FITTING OF THE ANGLE-DEPENDENT 
MAGNETORESISTANCE IN TL2201 

In this section we review how the Boltzmann trans- 
port equation can be used to fit the ADMR data. We 
begin with the simplest case whereby both lUc and r are 
isotropic before going on to discuss the more general case 
in which both parameters are anisotropic within the con- 
ducting plane. 
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FIG. 4: (Color online) Diagram describing the ADMR exper- 
imental technique, whereby the sample (schematically shown 
in blue) is rotated by an azimuthal angle (^iexp with respect to 
the laboratory frame {xp,yp, Zp), and then data are continu- 



ously taken as a function of polar angle 



In the actual 



experiment, the crystallographic {xc,yc) plane does not lie 
exactly in the laboratory frame due to a slight misalignment 
and the corresponding polar angle ^crys = (taken as the an- 
gle between the field direction and the normal to the plane on 



the sample) differs from 



Moreover, the azimuthal 



may change as a function of 6exp as explained in Appendix A. 
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where is the reciprocal space direction parallel to 
the magnetic field H, dip = ujcdt and lOq is considered 
isotropic. In order to use Eq. 5 to analyse the ADMR 
data, we follow Yamaji'^^ and define a vector fc^ as shown 
in Figure 5(b). The projection of the magnetic field on 
the azimuthal plane relative to the kx-Bxis (correspond- 
ing to the Cu-O-Cu bond direction) is labelled 0crys (see 
Figure 5 (a)). Each orbital plane is then defined by three 
parameters: the polar angle ^crys = ^, the azimuthal 
angle 0crys and fc^. The former two are determined dur- 
ing the experiment whereas the latter is an integration 
variable in the fitting procedure described below. The 
intersection of this plane with the FS gives the path of 
the quasiparticle in reciprocal space. This plane is given 
by the equation 



A. Isotropic T and ujc 

In the presence of a magnetic field a quasiparticle tra- 
verses the FS following the contours defined by the dis- 
persion. During this journey the quasiparticle will gain 
velocity from the electric field until it encounters a scat- 
tering event, after which it begins its journey again. As 
the angle of the field with respect to the crystal axes 
is adjusted, the quasiparticle will traverse different or- 
bits and the average velocity in the direction of the cur- 
rent can vary dramatically. This picture is formalized 
in the Chambers' tube integral, which is the solution to 
the Boltzmann transport equation in the relaxation-time 
approximation 

= ^ X/k (f ) ..(k,0) /\,(k,t)e*/^d, 

(4) 

where /k is the mean occupation of state k and r is as- 
sumed to be independent of k (or equivalently, isotropic 
in the azimuthal angle 0) . The Chambers formula can be 
used in a situation where both closed and open orbits are 
present. ■^■^'■^^ For our particular interest only closed orbits 
are involved (the FS is q-2D)2°'2^ and we are able to use 
the simpler Shockley-Chambers tube integral. Further- 
more it is easier to use cylindrical coordinates, in line 
with our description of kpln, 0). The interplane conduc- 
tivity is then given by 



kxSinO + kzC0s6 = \]s.b\ = fc°cos0. (6) 

This is a convenient notation because k^ can be 
uniquely described in terms of fc^ and the projection of 
the Fermi wave-vector onto the azimuthal plane as the 
quasiparticle traverses an orbit, k|,(0). In summary 



k, ^\kB\ = k°- 4((/))cos(0 - (/)crys)tan6». (7) 



We make the replacement e — > in Eq. 5 by ap- 
proximating df^ /de S{e — ep) for kT ^ ep- The 
periodicity of Vz{(t>) and Vz{(t> — (j)') in <j) and (f) — 4>' re- 
spectively is of some computational benefit. Taking the 
Fourier transform of v^'^'^'^^ and writing it as a Fourier 
sum gives 



Vz {(f)) = flo + flnCOS ncj) + 5„sin nt/), 



oo 

Vz{(j) - (j)') — co + ^ c„cosn((/) - (j)') + dnSmn{(j> — cj)'), 



(8) 

where an,bn,Cn and d„ are Fourier coefficients. Using 
a Laplace transform and after a few algebraic manipula- 
tions, the conductivity is finally given by'^^'^^ 
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FIG. 5: (Color online) Schematic representation of a single quasiparticle orbit with an oriented magnetic field. Panel (a) gives 
the projected view of the Fermi surface, defining the azimuthal angle 0crys, while panel (b) gives the definition of the parameters 
ks, k° and 6'crya. 
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where ujc = t^o Sind t — tq to emphasize that these 
parameters are isotropic. Eq. 9 is used to calculate the 
resistivity in the transverse direction pzz by taking the 
inverse of azz, which is correct to a good approximation 
due to the large anisotropy of the in-plane and interplane 
resistivity. Because parameters such as the effective mass 
m* are not well known, it is usual practice to simulate the 
relative change in magnetoresistivity Apzz/ PzzOi where 
Pzzo is the interplane resistivity at zero field, rather than 
Pzz directly. This normalization procedure means that 
the warping parameters in the kz direction can only be 
determined as ratios. In other words, the ADMR can be 
used to obtain values for fcgi/fei and kiQi/k2i but not 
^21,^61 or fcioi directly. 

The parameters we wish to determine therefore are 
^00, fco4, fcei/fei, ^101/^21 and ujqTo (the cyclotron fre- 
quency and the scattering time always appear as a prod- 
uct in the sum of Eq. 9 and thus behave as a single 
parameter). In a number of earlier studies on differ- 
ent T12201 crystals, in which all parameters were al- 
lowed to vary, a consistent set of FS parameters were 
obtained. ^^'^"'^^ This enables us to refine our parameter- 
ization and minimize the number of free parameters with- 
out losing confidence in their relative magnitudes. We fix 
koo for example by first obtaining the doping level p us- 
ing the universal phenomenological relation'^^ between p 
and the critical temperature Tc 



1 - 82.6(p- 0.16)^ 



(10) 



then adopting the simple hole-counting procedure, 



(7rA;2o)/(2V«)' = (l+p)/2. 



(11) 



Our next simplifying assumption is that t±{(j)) vanishes 
at eight symmetry points on the FS (see Figure 3) as 
expected from band structure calculations^'' and revealed 
by earlier ADMR measurements. For this to be the 
case, we require 



k2i 



fcioi 
k2i 



0, 



(12) 



which fixes fcioi/A;2i to whatever value fc6i/fc2i is given. 
Hence, only three parameters fco4, ^61/^21 and the prod- 
uct loqTq, are used to fit simultaneously five polar angle 
sweeps at different azimuthal angles (in other words, the 
data are treated as a single data set, not five separate 
curves). It is important to realize that these constraints 
could be relaxed without a significant effect on the other 
key parameters. 

The fitting procedure begins by evaluating Uzz for a 
given polar angle ^crys and azimuthal angle 0crys- The 
c-axis velocity is evaluated Vz — h'^^deik) / dkz as a func- 
tion of 4> for a given fc^, where e is defined by Eq. 3. The 
kz dependence is determined by Eq. 7 and substituted 
into Vz- For the given fc^, the Fourier transform is taken 
and the sum in Eq. 9 is evaluated. This is then inte- 
grated over k1 across the whole Brillouin zone and the 
result inverted to give Pzzi^crys, ^crys)- This is calculated 
for all ^crys and <?!>crys in a single data set. This process is 
repeated for different parameter values until a best fit is 
achieved using standard minimization procedures. 
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FIG. 6: (Color online) The solid lines are c-axis magnetoresistivity data Apzz = Pzz{H) — pzzo, normalized to the zero field 
resistivity pzzo, taken at diflerent azimuthal rotations (^crys = 8° (red), 0crys ~ 18° (orange), 0crys = 32° (green), i/icrys ~ 46° 
(blue), (jjcrys = 56° (violet) relative to the Cu-O-Cu bond direction) at three different temperatures (a) 4.2 K, (b) 14 K and (c) 
50 K. The azimuthal angles given here strictly apply only to Scxp ~ ±90 due to misalignment of the crystal with respect to the 
platform axes (see Appendix A for details). The black dashed lines are the best least-squares fits obtained assuming that ujcT 
(= uiotq) is independent of (j> and that the parameters fco4 and ^61/^21 are fixed to their values at T = 4.2K. Thus only ujqto is 
allowed to vary with temperature. 



The solid lines in Figure 6(a) are ADMR data taken 
at T = 4.2 K and iiqH = 45 Tesla, normalized to the 
zero field resistivity value Pzzo- Each color represents a 
different azimuthal angle at which the individual polar 
ADMR sweeps were taken. Despite the fact that ujqTq is 
less than 0.5 in this sample, the variations in the c-axis 
magnetoresistance are significant, both with azimuthal 
and polar angle, thus tightly constraining our parameter- 
ization. Note that these data were obtained on a different 
crystal to those reported in Refs.^"'^^ though the result- 
ing parameterization (fcoo = 0.729A~^, fco4 — -0.022A^^, 
^61/^21 — 0.7) is very similar. The best least-squares 
fits to Eq. 9 are shown as black dashed lines and appear 
quite adequate for the full range of azimuthal and polar 
angles studied. (Data at larger angles were not taken 
at this temperature in order to avoid the large torque 
forces that accompany a transition to the superconduct- 
ing state, which occurs here when Hc2{d) surpasses 45 
T.) 

Corresponding data and fits for T — 14K and T — 50K 
are shown in panels (b) and (c) respectively. For the 
fits at higher temperatures (where a larger angular range 
can be swept), all FS parameters are fixed to their 4.2 
K values and only the product ljqTq is allowed to vary. 
The fits rapidly deteriorate as the temperature is raised 
and are clearly no longer a reliable representation of the 
data. In fact, even if we allow fcoo,A;o4 and ^61/^21 to vary 
with temperature, the fits do not significantly improve. 
Furthermore, if fcoo is allowed to be a free parameter, 
the fitting procedure tends to minimize at values where 
the Fermi surface is larger that the first Brillouin zone. 



which is clearly unphysical. In response to this failing, 
we abandon our naive picture of isotropic ujc and r and 
proceed to incorporate anisotropy into the formalism. 

B. Isotropic r and anisotropic uJc 

To illustrate how significant anisotropy in lJc can be, 
we consider first the most elementary tight-binding de- 
scription of an isotropic square 2D lattice. The dispersion 
of such a system can be described by the equation 

£(k) — £0 = — 2t[cos(fc^a) + cos(fc2,a)], (13) 

where e — eo describes the quasiparticle dispersion taken 
relative to some reference (for example, the non-bonding 
energy So). Quasiparticles complete orbits with a fre- 
quency LOc that depends on the scalar product kp-vp via 
the expression 

..{c,.e)^eBcose^l^§-^. (14) 

Near the bottom of the band, the quasiparticle orbits in a 
magnetic field appear almost circular and vp is isotropic 
and nearly parallel to the crystal momentum k. As 
ep approaches the van Hove singularity (vHs) however, 
anisotropy in vp becomes significant^'^ and uJc develops 
four-fold anisotropy that essentially becomes infinite at 
the vHs. 
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Let us now turn to consider the analogous situation in 
T12201. According to recent ARPES experiments,^^ the 
FS can be fitted by a tight-binding dispersion relation 



e — So — ^(cosfca; + cosfcy) + t2{coskxCosky) 
+ ^-{cos2kx + cos2ky) + ^■{cos2kxCOsky + cos2kyCoskx) 
+t5{cos2kxCos2ky) , 

(15) 

with h ^ -0.725, t2 = 0.302, ts =0.0159, h = -0.0805 
and t5 = 0.0034(eV). 
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FIG. 7: (Color online) The parameters kp, vp and uJc for four 
different doping levels {p — 0.17, 0.22, 0.26 and 0.3) based on 
the dispersion relation given by Eq. 15 and assuming a rigid 
band shift. The black dashed line corresponds to the doping 
level of our Tc = 17K sample. 



In order to visualize the doping evolution of the FS 
parameters according to Eq. 15, we show in Figure 7 
the variation of kF(^), 'Vf{4') and ujc{4>) (left, center and 
right panels respectively) for different values of the chem- 
ical potential assuming a simple rigid band shift. The en- 
ergy contours have been centered on the X point of the 
Brillouin zone. Though the band structure may change 
with hole doping, this approximation scheme serves as 
a good illustration of how the anisotropy in ujc varies 
in a comparable way to that in vp. As expected, the 
anisotropy grows as the FS at (tt, 0) approaches the vHs, 
though in the doping range relevant to T12201, it never 
exceeds 30%. Interestingly, as the chemical potential is 
raised, the anisotropy of Wc changes sign, so that the cy- 
clotron frequency goes from being maximal to being min- 
imal along the zone diagonal, but retaining four-fold sym- 
metry throughout. We approximate this without making 
any assumptions as to the sign of lUc using the expression 



c^e(0)"' «wo"'(l+/3cos(40)). 



(16) 



The curves in the right-sided panel of Figure 7 corre- 
spond to a range of jS values from (3 ~ 0.3 (for p = 0.3) 
to /? « at optimal doping. 

The addition of this extra parameter causes only minor 
modifications to the fitting procedure. The conductivity 
is now replaced by the equation 
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xexp[ft,((/)) - h{(t) ~ 0')], 

(17) 

where /i(0) — — J ^j-jjyp^- Under isotropic circumstances 
h{(j)) = —(fi/uJoTo as in Eq. 5. However, in the case where 
LUc satisfies Eq. 16, this becomes 
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(18) 



We can now define two new periodic functions Pz {(f>) and 
Pz{4> — 4>') whereby 



(19) 



The Fourier transform of each function is given by Eq. 
8. The form of Eq. 9 is identical, only that is inter- 
preted as the average of ojc within the plane (see Eq. 16). 
The fitting procedure proceeds as described in Section 
IV A except our fitted parameters are now fco4, ^61/^21, P 
and luoTq. As before, koo = 0.729 A~\ ko4 = -0.022A-^, 
kei/k2i is fixed at low T (4.2 K) and only looTq and /3 
are allowed to vary as a function of temperature. Fig- 
ure 8 (a)-(c) shows the best least-squares fits of the same 
ADMR data under this new parameterization scheme. 
While the fits are closer to the real data than in the 
corresponding isotropic case, there is still a clear prob- 
lem with the higher temperature fits. If we choose to 
allow kQi/k2i to vary however, the fits become reason- 
able at all temperatures, as shown in Figure 8 (d)-(f)). 
The mathematical reason for this is that P has two com- 
peting roles: it appears in the exponent h{4)) and in the 
ratio Vz{(t>)/uJc{4')- In the latter, P plays a similar role to 
^61/^21, as can be seen with an expansion using elemen- 
tary trigonometric identities 



(sin2(/) -I- ■^^sin6(/) - 

k2i 



= 1-/3/2 



f3k 



61 



kei 
k2i 



/3/2 + 
fcioi ^ Pk^i 

k21 2fc21 



2fc21 

, PkiQi 



k2i 
sin20 



sinlO(?!))(l + /3cos4(/)) 



(20) 



2fc21 

sinlO(/i 



/3fci 



2fc21 



■sinl4(/). 



The ^61/^21, ^101/^21 and /3 terms can compensate 
each other as long as /3 remains small (that is, so long as 
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FIG. 8: (Color online) The raw data (solid curves with colors corresponding to ^crys as defined in Figure 6) and best least- 
squares fits (black dashed lines) for ADMR taken at three difi'erent temperatures T = 4.2K (panels (a) and (d)), 14K ((b) and 
(e)) and 50K ((c) and (f)). In panels (a), (b) and (c), the parameters fco4, fcgi/fei are fixed at their 4.2K values whilst t^oTo 
and /3 are allowed to vary. In panels (d), (e) and (f), kei/k2i is also allowed to vary with temperature. 



products such as /3fcioi/2fc2i are negligible). The only 
non-compensating contribution of (3 in this expansion 
is in the multiplication of the sin20 and sinl40 terms, 
though perturbations of the former would be noticeable 
first. In other words, the fitting procedure tends to keep 
the sum (3/2 + fc6i/^2i constant as a function of temper- 
ature and so the T-dependent changes are contained in 
the behavior of h((j)). We return to this point later in our 
discussion of Figure 10. 

The changes in kQi/k2i and P required to satisfactorily 
fit the data are significant (~ 20% change in A:6i/fc2i and 
a factor of 10 increase in /?) and if correct, would imply 
pronounced FS reconstruction with increasing tempera- 
ture. Between 4K and 50 K, one may expect the Fermi 
distribution to broaden by around 2% of the band width 
about the chemical potential. At p = 0.26 this is ap- 
proximately equivalent to a change in nominal doping 



of ±0.02 allowing a change in (3 of at most ±0.05, as 
is evident from Figure 7. This is significantly less than 
is required to quantitatively account for the evolution of 
the ADMR. To our knowledge, the FS restructuring re- 
quired to fit the present data has never been reported 
in cuprates and so justifying it would require some very 
subtle physical arguments. Indeed, photoemission stud- 
ies have reported insignificant changes as a function of 
temperature on overdoped compounds.**^ Moreover, in 
a recent doping dependence study, we found that the 
overall anisotropy decreases with increasing carrier con- 
centration, i.e. as one approaches the vHs, in marked 
contrast to the band structure picture discussed above 
and illustrated in Figure 7. In the following section there- 
fore, we turn to consider the effect of anisotropic scatter- 
ing, which not only allows the data to be fitted accurately 
but also avoids the physical and mathematical difficulties 
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we have encountered when considering anisotropy in Wc 
alone. 



C. Anisotropic r and anisotropic Wc 

In the most elementary description, the scattering Ufe- 
time r is the average time between collisions of an elec- 
tron travelling in a metal. In a FL picture however, this is 
taken to be the mean lifetime of an electron excitation, 
giving the decay time of a quasiparticle to its ground 
state near the chemical potential /y,. The rate of change 
of occupation of a state at k is related to the intrinsic 
transition rate between two arbitrary states k and k', 
weighted by the occupation of k and the lack of occupa- 
tion of state k'. 

The transition rate cannot be calculated without a pri- 
ori knowledge of the scattering processes that are present. 
In the limit of clastic scattering however, both k and k' 
arc on the same energy surface, and this function is sim- 
ply cos(^fefe'), where ^fefe' is the angle between kand k'. 
We then use the relaxation-time approximation, whereby 
Pkk' — P(0fcfc') is a function of (j)kk' only and this allows 
a natural definition for the scattering time r 



always appears in the product ujrjT in the sum of Eq. 26 it 
is clear that the procedure for incorporating anisotropic 
r will be identical to Section IV B where we introduced 
anisotropy in lOc- The simplest model would involve a 
four-fold anisotropy and in a similar fashion to Sande- 
man and Schofield^'^ or loffe and Millis^^ we write 

r(<^) =ro(l+acos(4(/))), (23) 
where Fq = I/tq. 




'x /x 

FIG. 9: The mean free path 1 for three different cases of 
the anisotropy parameter a for T12201 with p = 0.26: (a) 
the isotropic case (dashed) a — Q, {h) a = —0.1,-0.2,-0.3 
(solid, moving inwards along the zone diagonal) and (c) a = 
0.1,0.2,0.3 (solid, moving outwards along (tt, tt)). 



^(X j{l -cos(^fcfc'))P((/)fcfc')sin((/)fcfc')#fefc'- (21) 

The relaxation time approximation is often an excel- 
lent starting point for interpreting transport data and r 
is usually considered to be independent of momentum, 
both of the initial and final state. A more general theory 
however would allow Pkk' to be a function of the initial 
state k. In this instance, it may be assumed that the 
form of Eq. 21 stays very similar,^^ except that Pkk' is 
now k-dependent, and hence the replacement r — » r(k) 
needs to be made. In this case, it can be shown^^ that 



ffkr(k) 



ffk 

r(k) 



{gw- 3k'}IPkk'dk', 



(22) 



where gk = /k — /k and arc the probabilities 
of an electron occupying a state k in the presence of a 
field and in equilibrium respectively. r(k) is the scat- 
tering rate. Eq. 22 is general enough to include inelas- 
tic scattering mechanisms too (involving energy transfers 
< ksT for any given scattering event) and this is a direct 
consequence of the relaxation-time approximation. Such 
details would be normally be contained in ^k, and the 
Boltzmann equation would be very difficult to solve. In 
the relaxation-time approximation however, these details 
are deliberately ignored and all that is required for Eq. 22 
to hold is that the statistical ensemble of quasiparticles 
returns to equilibrium between collision events."*^ 

Under these circumstances, we are able to define an 
anisotropic scattering time that will enter all of our cal- 
culations of the conductivity. Since the scattering time 



Figure 9 illustrates the effect of this form of scatter- 
ing rate anisotropy on the mean free path 1 for the FS 
derived in Eq. 15 (assuming p = 0.26). If a = 0, the 
scattering rate is isotropic and I(k) simply follows the 
form of VF(k) (panel (a)). If a > (panel (b), F(k) is 
maximal in the direction parallel to the zone axes and 
competes with VF(k). If, on the other hand, a < 0, F(k) 
is maximal along the zone diagonals, the anisotropy in 
l(k) is enhanced in this direction. 

With this definition of F(0) , the conductivity is iden- 



tical to Eq. 17 but now we have h{(p) = J 



aJc(0)T(0) 



+ i(a + /3)sin4,^+ffsin80}, 

with the periodic functions Pz{<p) and Pz{(t> — 
redefined as 

UJc{<P) 

The conductivity is then given by 



and 



(24) 



b') now 



(25) 
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FIG. 10: (Color online) The quality of fit parameterised by 
the sum E(5^ as a function of f), which is given a specific value 
between -0.4 and 0.3. Conversely, the parameters uoqto, a and 
^61/^21 are free to vary. The parameter E(5^ has a broad 
minimum, indicating that a good fit can be achieved for a 
broad range of /3. Insets (a) and (b) show respectively the 
interdependence of ^61/^21 and a on /3 (see also Eq. 20 and 
24). The axes in (a) are shifted as described in the text. 



where = cjo/(l + Q:/3/2). Equipped with Eq. 26 we 
can now follow the procedure described above. The vari- 
able parameters are now fco4, fcgi/fc2i, a, /? and uj^to- As 
before we fix fcoo (= 0.729 1-^) and fco4 (= -0.022l-i) 
and fit the low temperature data with all other param- 
eters free to vary. In this instance however we are now 
over-parameterized since a and /3 can compensate each 
other to within a factor of ±a/3/4. As a consequence, 
one cannot accurately quote absolute values for each pa- 
rameter individually, but rather the sum a-\- [3 (see inset 
(b) of Figure 10). 

We can parameterize the quality of the fits by the sum 
of squared differences of the data from the fitted curve, 
which is denoted E(S^. As |/3| becomes large, the terms a 
and fcei/fc2i are no longer able to compensate and the fits 
decline in quality. However, as shown in Figure 10, there 
is also a broad flat region over which is minimized 
and one cannot pinpoint the exact value of [3. The axes 
in inset (a) are shifted so that it is apparent that the 
sum P/2 -\- kei/k2i is pinned to a value of about 0.72, a 
fact which continues to be true at higher temperatures 
no matter what one forces /3 to be. Similarly, the value of 
a -\- P is pinned to nearly zero at low temperature. From 



the low-T parameterization used previously to set the FS 
parameters, we can settle on a value of /3 « —0.1 ± 0.1 
which is comparable to that estimated from the ARPES- 
derived dispersion despite being opposite in sign. 

Panels (a), (b) and (c) of Figure 11 show the resulting 
fits to the new parameterization scheme, in which only 
a and looTq are allowed to vary with temperature, for 
T = 4.2K, 14K and 50K respectively. In contrast to pre- 
vious schemes, the quality of the fits are comparable at all 
temperatures, without the need for any variation in the 
other parameters. Hence, by introducing T-dependent 
anisotropy in the scattering rate, there is no longer any 
need to invoke FS reconstruction to account for the evo- 
lution of the ADMR data. We therefore conclude that 
this is the most elegant and physically realistic parame- 
terization scheme of all those considered here. 

The sign of a is found to be positive, indicating that 
scattering is weakest along the zone diagonals, as deter- 
mined previously by azimuthal ADMR measurements.^^ 
As the temperature is raised, a increases markedly. This 
implies that the anisotropy resides in the inelastic, rather 
than the clastic scattering channel. As with anisotropy 
in LUc, anisotropy in the scattering rate can be FS- 
derived, e.g. due to FS instabilities such as charge- 
density waves, spin-density waves or antiferromagnetic 
fiuctuations. Spin, charge, or indeed superconducting 
fluctuations all have specific momentum (and frequency) 
dependence that is peaked at (or in some cases, con- 
fined to) particular regions in k-space. Anisotropy in 
can also signify additional physics due, for example, 
to strong electron correlations near a Mott insulating 
state or anisotropic electron- impurity scattering. ^■^ The 
present analysis cannot of course reveal the microscopic 
mechanism of the anisotropic scattering itself, but can 
identify some important characteristics of the scattering 
mechanism, such as its magnitude or its symmetry. Sys- 
tematic measurements, e.g. as a function of doping and 
or pressure, would then allow a detailed comparison with 
the various theoretical proposals and thus help to reveal 
important hints as to its microscopic origin. 



V. CONCLUDING REMARKS 

In this paper we have set out a detailed formalism for 
incorporating in-plane anisotropy, both in the cyclotron 
frequency and in the transport lifetime, into the analysis 
of interlayer magnetoresistance of a q2D metal. The fo- 
cus of the present paper has been to illustrate the need 
to introduce an anisotropic scattering rate in order to 
explain the evolution of the ADMR data in overdoped 
superconducting T12201 within a Boltzmann framework. 
An anisotropic cyclotron frequency Lud^j)) can fit the data, 
but only if we allow the parameters describing the Fermi 
surface itself to change as a function of temperature. 
Given the absence of evidence for such reconstruction, 
this hypothesis seems unlikely. If, on the other hand, 
an anisotropic scattering time is introduced, all the FS 



11 




0.0 



-80 



-40 



40 



0.0 



-40 





^exp 



40 



0.0 



-40 





'^exp 



40 



FIG. 11: (Color online) Raw ADMR data (solid curves with colors corresponding to </>crys as defined in Figure 6) plotted with 
the best least-squares fits (black dashed lines) for three different temperatures (a) 4.2 K, (b) 14 K, (c) 50 K. Here both ujc and 
T are considered anisotropic in the azimuthal angle (f>. The parameters ^04,^61/^21 and /3 are fixed to their values determined 
at T = 4.2K whilst luoto and a are allowed to vary as a function of temperature. 



parameters can remain constant and only wqTo and the 
anisotropy in T{4>) are adjusted. Such a simple param- 
eterization is both elegant and experimentally accurate 
and we therefore believe it to be the most likely explana- 
tion of the observed ADMR data. 

A cautionary note is perhaps appropriate here. In the 
preceding calculations we have assumed the relaxation- 
time approximation and so the microscopic relaxation 
dynamics have been ignored. The concept of anisotropic 
scattering remains valid as long as t((/)) is interpreted as 
the exponential decay of the distribution between scat- 
tering events'*^''*'*. However, we cannot rule out more 
exotic relaxation dynamics that depend on the presence 
of the magnetic field and hence may be manifested differ- 
ently if probed by a different means (for example, pho- 
toemission). We have discovered that such exotic dy- 
namics do not need to be invoked to explain our ADMR 
data and conclude that its evolution with temperature, 
when viewed from a Boltzmann framework using the 
relaxation-time approximation, is best explained by a 
scattering rate with a temperature dependent anisotropy. 

At high doping levels, lifetime separation in cuprates 
is less apparent,^* leading some researchers to consider 
the problem from this perspective. This route has the 
added advantage of allowing the limits of the conven- 
tional Boltzmann transport theory to be explored as one 
moves across the phase diagram towards to more ex- 
otic and potentially non-FL ground state on the under- 
doped side. The key message here is that by generaliz- 
ing the theory to include an anisotropic scattering rate 
r(k) = T^^(k) one can continue to apply the Boltzmann 
approach and successfully account not only for the evo- 
lution of the ADMR with temperature, but also the dis- 



tinct T-dependencies of pab and cotGn found in over- 
doped T12201.^^ Furthermore, initial measurements of 
the doping dependence of a in T12201 suggest a signif- 
icant increase in anisotropy in r^^(k) as one move to- 
wards optimal doping, consistent with the observed in- 
crease in lifetime separation (as manifest in the temper- 
ature dependence of the Hall coefficient) with decreas- 
ing doping. ^'■^^'^^ The introduction of such anisotropy has 
proven a fruitful model to understand the normal state 
of high-temperature superconductors, though 
clearly more work is needed to parameterize r~^(k) fully 
and to identify the origin of the anisotropy. 

Finally, although the focus of this paper has been a sys- 
tem with body-centered-tetragonal symmetry, the anal- 
ysis could very easily be generalized to layered systems 
of other crystallographic symmetries as already pointed 
out in Rcf.^"^ and perhaps also one-dimensional sys- 
tems with anisotropic scattering. ADMR experiments 
on BEDT-TTF based organic superconductors have al- 
ready been performed at low temperature and explained 
in a Boltzmann framework without the need to invoke 
an anisotropic scattering rate.'^^ A full azimuthal and 
temperature dependence on other salts may however re- 
quire the introduction of such a parameterization. Simi- 
larly the same ideas may also apply to layered charge- 
density-wave compounds such as the rare-earth tritel- 
lurides RTe3.^° The Boltzmann equation, though simple 
in its assumptions, thus remains a powerful paradigm 
whose explanatory power is still to be explored. ADMR 
is an ideal probe for just such an exploration. 
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APPENDIX A: ACCOUNTING FOR SAMPLE 
MISALIGNMENT IN THE ADMR FITTING 
PROCEDURE 

In this section we illustrate how sample misalignment 
can be accounted for. The effect of sample misalignment 
on ADMR has been considered by several authors before 
this study, in particular Goddard^-'^ and Abdel-Jawad .^"^ 
Before we begin, let us define two frames of reference, 
that of the laboratory Xp,yp, Zp, in which the field is par- 
allel to the Zp direction, and that of the crystal Xc, Uc, -Zc- 
The Up axis is taken as the axis of polar rotation and 
Xp is perpendicular to this. The normal to the crystal- 
lographic plane shall be defined as Zc, and the in-plane 
directions Xc and Pc shall be taken to be parallel and per- 
pendicular to the copper-oxide bonds respectively. The 
angle between the field direction and the crystallographic 
normal Zc gives the crystallographic polar angle 9. The 
projection of the field onto the Xc — Vc plane gives the 
crystallographic azimuthal angle (j), taken from the Xc 
axis. 

There are two important differences between this study 
and that of Goddard et al..^^ Firstly, instead of correct- 
ing experimental ^exp and ^exp for misalignment to find 
the appropriate crystallographic 6 and (j), we fit the ex- 
perimental data by including the misalignment in the fit- 
ting procedure. Secondly, in the analysis of Goddard, '"'^ 
the crystallographic axes Xc and Uc can fall anywhere in 
the plane of the crystal and do not have assigned direc- 
tion with respect to the crystal bonds. This gives the 
misalignment one less parameter, and one of the crystal 
axes can always fall somewhere in the Xp — yp plane of 
the laboratory frame. In the present analysis this is not 
the case and so three rotations need to be included in the 
fitting procedure in order to account for every possible 
misalignment: 

• Beginning with the sample aligned with the labo- 
ratory frame, there is a misalignment in azimuthal 
angle denoted by $asym, rotated about Zp as shown 
in Figure 12 (a). 

• A rotation about the yp axis denoted by QX^yai ^ 



shown in Figure 12 (b). 

• A rotation about the Xp ctXlS denoted by O^ym as 
shown in Figure 12 (c). 

These transformations are elegantly described alge- 
braically. We follow the notation whereby a rotation 
Rf){6)a is a rotation of a vector a by angle 6 about an axis 
[3. In particular the laboratory axis z'^ is tranformed rela- 
tive to the crystallographic Zc{Q^ 0) (before any azimuthal 
or polar rotation) axis to 

4 = R-^i^t.yJRyi,iQlyMO, 0), (Al) 

due to the misalignment of the sample. In an ADMR 

experiment, the sample is then rotated about the labora- 
tory Zp azimuthally by an angle (f>cxp, and then rotated 
about the polar angle Ocxp- The position of the 

crystallographic Zc(6'cxp, '/'cxp) after these azimuthal and 
polar rotations axis is given by 

-^-C (^Cxp 7 4^CKp 

With reference to Figure 13 it is elementary to see that 
the projection of the field parallel to Zp on the crys- 
tallographic Zc(^cxp, (^oxp) will give ^crys, which should 
be used to calculate the value of the magnetoresis- 
tance in the analysis. Similarly, the projection on the 
2:c(6'exp,</'exp) - 2/c(6'crys,</'exp) plane will yield 0crys- Al- 
gebraically we have, 

COs(^crys) = Zp ■ Zc{deKp, (/"exp), (A3) 

and 



The asymmetries for the sample considered in the present 
paper were estimated to be <I>asym = 8°, 0|gy„j = 0°, 
®asyin = 3°. The parameter $asym was approximately 
equal to values estimated from diffractometry performed 
after the ADMR experiment. 
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FIG. 12: (Color online) The reference frame of the laboratory Xp,yp,Zp and that of the crystal Xc,yc,Zc where Xc is parallel 
to the copper-oxide bonds, (a) The first misalignment considered is that where the experimental Zp and crystal Zc directions 
coincide but the axes on the azimuthal plane are offset by an amount "^asym. The second misalignment considered is that where 
the experimental and crystal z directions do not coincide. Two rotations are responsible, one about the platform j/p axis by an 
angle Basym shown in (b), and another about x'p by an angle Oasym shown in (c). This gives completely general description of 
the crystal misalignment with respect to the platform axes. 
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FIG. 13: (Color online) The polar Scrys = and azimuthal 
4>crys angles that enter Eq. 17 depend only on the projection 
of Zp on the crystal axes. This projection will change as a 
function of ^cxp and ^cxp and so must be calculated at each 
point. 
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